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In this work, a cosmological model inspired by string theory with Gauss-Bonnet term 
coupled to the fermionic field is taken into consideration. The self-interaction potential 
is considered as a combination of the scalar and pseudo-scalar invariants. Here the cos- 
mological contribution of the coupling of Gauss-Bonnet term with a non-Dirac fermionic 
field-characterized by an interaction term L 2 d q— is investigated. It is observed that the 
new type of coupling plays a significant role in the accelerating behavior of the universe. 
Specifically, in addition to the late time acceleration for the universe, L 2 DG produces an 
early decelerating behavior. The behavior of the equation-of-state parameter (w) is such 
that it guarantees the stability of the theory. 
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1. Introduction 

Nojiri et. proposed the Gauss-Bonnet dark energy model, inspired by string/M- 
theory where standard gravity with scalar contains additional scalar-dependent cou- 
pling with Gauss-Bonnet invariant. Their study indicates that current acceleration 
may be significantly influenced by string effects. Also, the current acceleration of the 
Universe may be caused by a mixture of scalar phantom. Moreover, they concluded 
that scalar-Gauss-Bonnet coupling acts against the occurrence of Big Rip in phan- 
tom cosmology. Ribas et.al^l, in 2005, investigated whether fermionic field, with a 
self-interacting potential that depends on scalar and pseudo-scalar invariants, could 
be responsible for accelerated periods during the evolution of the Universe, where 
a matter field would answer for the post inflation decelerated period. They have 
shown that the fermionic field behaves like an inflation field for the early Universe 
and later on, as a dark energy field, whereas the matter field was created by an 
irreversible process connected with a non-equilibrium pressure. 

It is expected that the combination of fermionic field and Gauss-Bonnet term 
will produce significant effects in cosmological evolution. Therefore, we propose 
a cosmological model with fermionic field and Gauss-Bonnet term. The coupling 
of Gauss Bonnet term to the fermionic field is introduced in order to produce a 
late time acceleration for the universe. We explore whether the combined effect 
can produce an early inflation as well as late time acceleration. Taking inspiration 
from string/M-theory in which additional higher order terms are required, a term 
proportional to the square of the fermionic Lagrangian density which we call as an 
interaction term, is added to the total Lagrangian. In this model, the additional 
coupling is suppressed by some powers of the string scale. So we produce a long 
term effect without disturbing the overall stability. It may be noted that the fermion 
condensate expected from Super Symmetry theory also calls for higher order terms. 

2. String Theory Inspired Cosmological Model with Fermionic 
Field and Gauss-Bonnet Term 

The action for this model reads 




where L g = R/2, with R denoting the curvature scalar, L g is the gravitational 
Einstein Lagrangian density L rn is the Lagrangian of the matter field and aiand a 2 
are coupling constants. Finally, the augmented Dirac Lagrangian density Lug for 
a fermionic mass m is given by 

Log = \ [Wr^* - (£> M ¥) r"*] - m - VL GB (2) 

In Eq. ([2]), V = V(^f, ^f) describes the potential density of self-interaction between 
fermions, Lqb is the Gauss-Bonnet Lagrangian density^, given by Lqb — R 2 — 
AR^RV + R f , l/p(7 R^ f " J and * = 4^7° denotes the adjoint spinor field. Moreover, 
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the connection between general relativity and Dirac equation is done via the tetrad 
formalism and the components of the tetrad play the role of gravitational degrees 
of freedom. That is, — e%j a are the generalized Dirac-Pauli matrices, where 
denote the tetrad or vierbein. The covariant derivatives ^ in Eq. ((2|) are given by, 



= + ttfi M (3) 

with the spin connection, 



1 

^ypvi^ ps ~ c 6 

where T^ s is the Christoffel symbol 



n, = --g pa [K s ~e p b d,e b s }^r s (4) 



2.1. Basic Equations 

Based upon the cosmological principle, which states that our universe is homogenous 
and isotropic, we use the Robertson- Walker metric to describe our Universe: 

ds 2 = dt 2 - a{t) 2 (dx 2 + dy 2 + dz 2 ) (5) 

Through Euler-Lagrange equations, from Eqs. ((!]) and @, we can obtain the 
equations of motion for the spinor field as 

(oi + 2a 2 L DG ) dL DG /(ty = (6) 

which requires that either (ai + IhiLdg) = or dLrja/dTp = 0. 

The first case, i.e. (ai + 2o,iLt)g) = , gives Log = — ai/2a 2 . Being forced to 
have Ljjg equal to a constant is an unphysical and artificial constraint. Hence we 
rule out that option. 

Now choosing the second case, we get the Dirac equations for the spinor field as 
dL dg /dip = or 

iT»D^-m*-L GB ^==0 (7) 

The equation of motion is exactly the same as in the case where no term is 
present. However we see that the L 2 DG term modifies the source terms for gravity in 
a significant way. We also note that any series of arbitrary powers L DG present in 
the original Lagrangian would lead to SLdg/SiP = for the fermionic equation of 
motion, because the other equation arising would be unphysical. So our analysis can 
be thought of as treating the case of a general form for the fermionic contribution 
to the stress tensor, where we retain terms only upto the next to leading order 
contribution in an expansion of the stress tensor in the derivatives of ip. 

Similarly, we can obtain the Dirac equations for the adjoint coupled to the 
gravitational field. 



October 20, 2009 8:14 WSPC/INSTRUCTION FILE Rakhi'final 



4 Rakhi et al. 



Also the variation of the action (Eq. (JTJ) ) with respect to the tetrad gives the 
Einstein's field equations as 



(8) 



where T^ v = Tf + T^ 1 with Tf being the energy- momentum tensor of the 
fermionic field and T^ , that of the matt er fi eld. The symmetric form of the energy- 
momentum tensor of the fermionic field which follows from Eq. |T]) gives 

Tf = - {^VD 1 '^ + ^T U D^^ - D"*^* - D^T l '^} + -L DG {^T^D u ^^T u D'"^- 



D V ^T^ - D^VV^} - g^L DG - g^L 2 DG 



(9) 



where the Dirac matrices have the form T = 7 , T l — ^75-7*, T 5 = 7 s and the spin 
connection components are 



1 



i 



= 0, a = - a(t) 7 l 7 



(10) 



From Eq. ([9]), together with Eqs. JTJ), © an d ©, we can obtain the non- vanishing 
components of the energy-momentum tensor for the fermionic field as 

(T/)° = ax (m(^) + VL GB )~a 2 (m (¥ *> + VL GB ) 2 +a 2 fm(W*) + ^^L GB + ^^L GB 

(11) 



{T f )\ = (T f )l = (7»1 = ai L GB 



V 



dV dV^> 
~2 ~ d¥"2 



-a 2 L 



2 

GB 



V 



^ dV dV^ 



2# d^ 2 
(12) 

The interaction between the fermionic constituents are modeled through the pres- 
ence of a non-equilibrium pressure term(-n7) in the source's energy-momentum ten- 
sor. Therefore, for an isotropic, homogenous Universe, the components of the energy- 
momentum tensor of the fermionic field can be identified as 



{ T f)u = diag (pf,~Pf, ~Pf, -pf) ■ 



(13) 



One can obtain from Eqs. (fTTj) , (fi"2|) and (fl"3j) . the energy density and pressure of 
the fermionic part: 



p f = a 1 (m(^)+VL GB )-a2(m(^)+VL GB ) 2 +a 2 (m(^) + ^^L GB - 



(14) 



'GB J 



Pf = aiL GB 



V 



* dV dVf 



2 rf* 2 
The self-interaction potential takes the form, 



a 2 L 



2 

GB 



V 



* dV dV^ 
~2 M _ rf*"2 



V = A 



(15) 



(16) 
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where A is the coupling constant and n is a constant exponent. We consider V as a 
combination of the scalar and pseudo-scalar invariants. That is, /3i = /?2 = 1. 

Combining Eqs. (| 14[) . (fl~5f and (|T6|) we solve for the energy density and pressure 
of the fermionic field as 

p f = a\ (m + UL GB ) - a 2 (m (**) + VLgb) + a 2 (m + 2nVL GB ) 

(17) 

P/ = ai (2n - 1) V Lqb + a 2 [(2n - 1) Fics] 2 (18) 

respectively. We infer from Eq. fTSf that the pressure is negative only for n < 1/2 
and this could represent either the inflaton or the dark energy. 
The Friedman equations give 

■m 2 = P = p f + Pm (19) 

2- = -p- w = -H 2 -p- w (20) 

Qj O 

or 

2^ + i? 2 = -p/-p m -w (21) 

Using Eqs. I|19p and (|2ip. the acceleration equation becomes 

» = 2ff 2 - p f - p m - Pf - Pm -zu 
a 2 

The conservation law for the energy density of the fermionic field and matter field 
can be written as 

Pf+3H{ Pf +p f )=0 (23) 
p m +3H { Pm + Pm + w) = (24) 

where H = ^^-is the Hubble's parameter. 

Considering the fermionic field to be a function of time alone, the Dirac equations 
([5|) becomes 

*+^i/* + iTO7°^ + i7° J L Gi 3^ =0 (25) 
2 dvp 

Similarly, the Dirac equation for the adjoint spinor field coupled to the gravitational 
field becomes 

2 ' <m 



2.2. Field Equations 

The system of field equations that help us to find the cosmological solutions is: 
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2.2.1. The Acceleration Equation 

The acceleration equation for the cosmological model under consideration can be 
obtained from equations (jTTJ) , (JTSJ) and (|2"2")l . Also we consider the pressure of the 
matter field as p m — w m pm with < w m < 1, which is a barotropic equation of 
state. We again consider a vanishing energy density of the matter field (p m ) and 
a vanishing non-equilibrium pressure (w)a,t t=0. Incorporating all these ideas, the 
acceleration equation for the present model becomes 



a 
a 



c±d 



where 



and 



<■ = I °: °L (1 + 7.2 ai H 2 V - 9.6a 2 H 2 Vm (##)) 



2 a 2 H 4 V 2 



(27) 



(28) 



d = l\j^§syT + 7-2 ai HW - 9.6a 2 H 2 Vm (W)) 2 - ^]y 2 (2# 2 - a x m (**) - p m - w m p m - w) 

(29) 



2.2.2. Evolution equation for the energy density of matter field 
The evolution equation for p m is given by 

p m +3H(p m + w m p m + w) = 



(30) 



2.2.3. Evolution equation for the non- equilibrium pressure 

The linearised form of the evolution equation for non-equilibrium pressur d^^ reads 



t m +zu + 3rjH = 



(31) 



2.2.4. The Dirac Equation 

In terms of the spinor components,^ = *2i *3i *4) T i the Dirac Eq. (123)) can 



be written as 













(Vx \ 


d 


*2 




*2 


+im 


*2 


dt 


*3 


*3 

uj 


-*3 
\-*4/ 


-48 


i±H 2 

a 




+ *J* 2 - * 


1*3 - * 



l^F 2 (V^l + * 2 *2 - *3*3 - *4*4) 

(32) 



*2 

-*3 

V-* 4 y 



dV 



/*3 \ 
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2.2.5. Equation -of -State Parameter 

The equation-of -state parameter (w) is given by 

w = V - (33) 
P 

This should be thought ol as a phenomenological relation reflecting the current 
amount of pressure and energy density in the dark energy. 

2.3. Cosmological Solutions 

In order to obtain numerical solutions of the coupled system of Eqs. (f2"T)) - (|3"3"|) , we 
have to specify the parameters A, /?2, n, m, a, w m and a\. These parameters take 
the following values: 

A = 0.1, (3\ = 02 = = 0.3, m = 0.01, a = l,w — 1 and a± = 1. The initial 
conditions, we have chosen for t=0 (by adjusting clocks) are: a (0) = l,$i (0) = 
O.li, * 2 (0) = 1, * 3 (0) = 0.3, * 4 (0) = i, p m (0) =0,w (0) = 0. The last two initial 
conditions correspond to a vanishing energy density of the matter field and a van- 
ishing non-equilibrium pressure at t = 0. The conditions chosen here characterize 
qualitatively an initial proportion between the constituents in the corresponding 
era. 

Graphs are plotted for various values of 02- The effect of L 2 DG is now studied 
first by choosing 02 = 1 [Figs. 1-5]. 

2.3.1. Acceleration 

In Fig 1, we plot the accleration as the root of {57]), corresponding to the upper + 
sign, whereas in Fig 2, the acceleration corresponding to the lower - sign is plotted. 




A3 OA CLB CLS 1 IS 1A 

rnne 



Fig. 1. 

From the following Fig 2, it is clear that the negative root has no interesting 
features. Hence from now onwards we concentrate only on the root corresponds to 
the + sign. 
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2.3.2. Energy density 

In Figs. 3 and 4, we plot the evolution of energy densities of fermionic field and 
matter field respectively. 
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Fig. 3. Energy density of fermionic field P/Vs time 




Fig. 4. Energy density of matter field p m Vs time 
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2.3.3. Equation- of- State Parameter 

The evolution of cquation-of-state parameter is plotted below: 

o 

-0.05 

■0.1 
-0.15 

-0.2 
-0.25 

-0.3 
-0.35 

Fig. 5. Equation of State Parameter of the universe Vs time 

The effect of L 2 DG is studied below, by putting a 2 = 5, 10, 50,100 successively 
[Figs.6-9]. 




Fig. 6. Acceleration Vs time for various <i2 
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Fig. 7. Energy density of the fermionic field P/ Vs time for various ai 
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O 1 2 3 4 



Fig. 8. Energy density of the matter field p m Vs time for various d2 
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Fig. 9. Equation-of-state parameter of the Universe Vs time for various 02 



2.4. Discussions 
2.4.1. Acceleration 

From the above obtained graphs, it is seen that the coupling of Gauss-Bonnet term 
with non-Dirac fermionic field produces an early inflation preceded by a short term 
deceleration. At later times, there is an increase in the acceleration effect. This 
late time acceleration is considered to be the effect of Gauss-Bonnet term. Also as 
the value of the coefficient of L 2 DG: 0,2 increases, the initial acceleration decreases 
and hte late time acceleration becomes more delayed. This shows that L 2 DG has 
a prominent role in the present model. It is possible to get the purely non- Dirac 
fermionic field, by simply putting a\ = 0. 
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2.4.2. Energy density / Equation of state parameter 

It is evident from the plots that the energy density of the matter field evolves more 
rapidly for smaller values of &2- As a2 increases, the late time evolution obtained 
for p m becomes more and more delayed. Also it is observed that the decay of the 
energy density of the fermionic field comes earlier for smaller values of a 2 , causing 
larger accelerated period. As the value of a-i increases, the decay of pf starts at a 
later time. 

These effects can be explained on the basis of equation-of-state parameter. The 
equation-of-state parameter is connected directly to the evolution of the energy 
density, and thus to the expansion of the universe. It is seen that the equation-of- 
state parameter is actually increasing during the initial stages of evolution. However, 
at later times, the equation-of-state parameter decreases and becomes negative. 
Finally, at very late times, it increases and saturates to a negative value. It is 
observed that w > — 1. This guarantees the stability of the theory and the saturation 
towards w > — 1 is considered to be resulting from the coupling of Gauss-Bonnet 
term with the non-Dirac fermionic field. 



2.5. Final Remarks and Conclusions 

In this work we have considered a cosmological model inspired by string theory in 
which Gauss-Bonnet term is coupled to the non-Dirac fermionic field. Here the self- 
interaction potential is considered as a combination of the scalar and pseudo-scalar 
invariants. 

We have investigated the possibility whether the coupling of Gauss-Bonnet 
term with a non-Dirac fermionic field-characterized by an interaction term L 2 DG ~- 
contributes significantly to the cosmological evolution. We have observed that the 
new type of coupling plays an important role in the accelerating behavior (both 
early and late time) of the universe. It is seen that the behavior of the equation- 
of-state parameter (w) is such that w > — 1, which guarantees the stability of the 
theory. 

Hence we conclude that the coupling of Gauss-Bonnet term with the non-Dirac 
fermionic field has a significant effect on the cosmological evolution. An interesting 
feature about this model is that this can generate standard acceleration similar to 
Dirac as well as a non-Dirac behavior by simply adjusting the magnitude of the 
coupling constant a.2- 
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